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We demonstrate, with a concrete example, how the combinatorial approach to a general system of 
particles, which was introduced in detail in the earlier paper [1| , works and where it enters to provide 
a genuine extension of results obtainable by more traditional methods of statistical mechanics. To 
this end, an effort is made to study cluster properties of the one-dimensional lattice gas with nearest 
neighbor interactions. Three cases: the infinite temperature limit, the range of finite temperatures, 
and the zero temperature limit are discussed separately, yielding some new results and providing 
alternative proofs of known results. In particular, the closed-form expression for the grand partition 
function in the zero temperature limit is obtained, which results in the non-analytic behavior of the 
grand potential, in accordance with the Yang-Lee theory. 
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I. INTRODUCTION 

In our recent paper [l| (hereafter referred to as the 
paper I), we have started a new line of theoretical re- 
search on imperfect gases and interacting fluids. We 
have exploited concepts of enumerative combinatorics 
to deal with equilibrium systems in the infinite volume 
limit. The approach seems to provide precise mathe- 
matical techniques relevant to study phase transitions. 
In particular, we have shown that the so-called perfect 
gas of clusters model underlying various cluster /droplet 
theories of phase transitions (see eg. Q and references 
therein) naturally emerges from our approach. In the 
model, the basic idea is that an imperfect fluid, which 
is made up of interacting particles, can be considered as 
an ideal gas of clusters at thermodynamic and chemical 
equilibrium. There is no potential energy of interaction 
between clusters and the clusters do not compete with 
each other for volume. The main conclusion drawn from 
the earlier paper was that the grand potential (the Lan- 
dau free energy) of such a clustered system of interacting 
particles may be considered as the exponential generat- 
ing function for the number of internal states (thermo- 
dynamic probability) of these clusters. 

In this paper we will use the combinatorial approach 
described in the paper I to analyze properties of the 
one dimensional lattice-gas with nearest neighbor inter- 
actions. The aim of this paper is twofold. First, we want 
to show on a specific case, how our approach to interact- 
ing fluids works, and where it enters to provide a genuine 
extension of results obtainable by more traditional meth- 
ods. Second, we want to demonstrate with a first con- 
crete example that the approach may provide important 
insights into microscopic mechanisms which lead to the 
occurrence of phase transitions. Although phase transi- 
tions in the traditional sense are proved not to exist in 
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the one-dimensional lattice gas considered in this paper, 
its behavior at T = OK seems to have some direct bear- 
ing to the problem. To show this, we will focus on the 
relation between thermodynamics and cluster properties 
of the lattice gas. The general question of such a relation 
has been raised by various theories of phase transitions, 
in which geometric interpretation of the liquid-gas transi- 
tion was shown to consist in the sudden formation of the 
macroscopic cluster. In the following, with the example 
of the one-dimensional lattice gas, we will show that our 
approach is well suited to handle such problems. 

The outline of the paper is as follows. In section [TTJ 
we briefly review results of the paper I. Section IHII is 
devoted to a detailed description of the one-dimensional 
lattice gas model and its equivalence with the chain of 
Ising spins in the external magnetic field. In section IIV1 
methods described in Sec. [TTJ are used to study cluster 
properties of the lattice gas. The paper is finalized in 
section [V] with some concluding remarks 



II. MICROSCOPIC MEANING OF GRAND 
POTENTIAL 



A. Combinatorial approach to a general system of 
particles 

In the paper I, we have considered a general system 
of interacting particles. The thermodynamic state of the 
system was given by the temperature, T, and the chemi- 
cal potential per molecule, \i. Assuming that the classical 
treatment is adequate, we have used the grand canoni- 
cal ensemble to describe the open system in the infinite 
volume limit, V — > oo. We have shown that the grand 
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partition function, 
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where 8 = (fc^T) -1 , z = e^, Z(B,N) is the canonical 
partition function of the considered system with N par- 
ticles, Z(/3,0) = 1 represents the so-called vacuum state, 
and g(E, N) stands for the density of states, can be writ- 
ten as 
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where the coefficients w n (8) are given by the par- 
tial derivatives of the grand thermodynamic potential, 
$03, z), i.e. 



w n {8) = -pct> n {(3) = -8 



dz r - 
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while B N>k ({w n }) = B Nik (wi, w 2 , . . . , w N - k +x), and 
B N ({w n }) = ELi b nAM) (with B ({w n }) := 0) 
represent the so-called incomplete and compete Bell poly- 
nomials, respectively. 

The partial Bell polynomials in Eq. ^ are defined as 
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where the summation takes place over all non-negative 
integers c„ > 0, such that the constraints J2 n c„ = k and 
nCn = N h°ld. Given combinatorial meaning of Bell 
polynomials, we have shown in the paper I that in the 
case when the coefficients w n (B) are all non-negative, i.e. 



V n >i w n (B) > 0, 



(9) 



the equality of the two series, Eqs. ([2]) and ([6]), give rise 
to the following formula: 

1 
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( 10 ) 

where f(k,E) is the probability that the system hav- 
ing energy E consists of k independent clusters. Also, 
we have shown that Eq. (| 10[) describes probability that 



a general system of N interacting particles at tempera- 
ture T, regardless of its energy, consist of k clusters. It 
has been argued that the formula provides a novel and 
strictly microscopic understanding of the grand thermo- 
dynamic potential, as the exponential generating func- 
tion for the numbers of internal states, w n (j3), of n- 
clusters, i.e. 



*GM = -~lnS(/?,*) 
w n (8) 
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Finally, in the earlier paper, we have noted that com- 
parison of Eqs. ([I]) and (O leads to a neat expression for 
the canonical partition function: 



Z(B,N) = ±^B N ({w n (j3)}) 

N 
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The expression is valid regardless of whether the condi- 
tions for the perfect gas of clusters model, Eq. ([3]), are 
satisfied. 



B. Perfect gas of clusters model 

The microscopic meaning of the grand potential de- 
scribed in the paper I directly relates to the perfect gas 
of clusters model (which is known from classical theory 
of simple fluids Q). In the model, one deals with a col- 
lection of non-interacting clusters, without knowing how 
one can build them as disjoint sets of interacting parti- 
cles. In the infinite volume limit, V — > 00, the pressure, 
the density, and the cluster size distribution (that is the 
mean number of clusters of size n) are respectively given 
by: 



and 
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N n (8,V) = Z n (8-,V)z n , 
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(17) 



where Z n (8,V) stands for the partition function which 
characterizes clusters of size n. 

The mentioned relation between our results and the 
classical gas of clusters model [22[ can be seen by first 
rewriting Eq. (|15j) in the following equivalent form: 
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and then comparing Eq. (| 18[) with the expression for the 
grand thermodynamic potential, Eq. (fT2|) . which raises 
the microscopic meaning of z). From the compar- 
ison it becomes evident that the conditions given by 
Eq. (|9]) make our approach to a general system of inter- 
acting particles equivalent to the perfect gas of clusters 
model with 



lim Z n (p,V) 
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III. 



ONE-DIMENSIONAL LATTICE GAS: 
SURVEY OF KNOWN RESULTS 



Let us consider a one-dimensional periodic lattice that 
consists of V sites (in the model, V represents volume) 
and a collection of N particles. The particles occupy sites 
of the lattice with the restriction that not more than one 
particle can occupy a given lattice site and only parti- 
cles on nearest neighbor sites interact. If we introduce 
variables Uj for each lattice site i, such that u\ = +1 if 
the site is occupied and <jj = otherwise, then the total 
energy for a given configuration of particles, {cr,}, is 



e g(Wi}) = -eyio-jO-i+i, 
»=i 



(20) 



where — e is the interaction energy between two neigh- 
boring particles, the periodicity of the lattice is imposed 
by assuming that ay+i — &i, and 



(21) 



In order to apply the general method described in 
Sec.|TT]to investigate the one-dimensional lattice gas, the 
first task is to find its grand partition function in the 
infinite volume limit, V — > oo, i.e. 



HgCS, V, z) = 1 + z N Z G (f3, V, N) 
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where Z G (f3,V, N) in Eq. (|22|) is the canonical parti- 
tion function for this gas and the starred configurations, 
{cr,:}*, in the second sum of Eq. (|2"3")l are those for which 
the condition given by Eq. (f2~Tj) holds. To carry on we 
exploit the mathematical equivalence between the grand 
partition function for the lattice gas and the canonical 
partition function for the Ising model in the external 
magnetic field [3]. 



The one dimensional Ising model consists of a chain of 
V spins, Si — ±1, with the interaction energy in a given 
configuration, {st}, given by 



El({Si}) = - J ^S,;S i+ i - i?y^Sj, 



(25) 



where J is the coupling constant between nearest neigh- 
bors, H is the external magnetic field, and sy+i = si. 
The canonical partition function for the model can be 
written as 



Z!(p,V,H) = ^ , 
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To show that Eqs. (|2"4")l and (|27p are in fact equiva- 
lent one has merely to note that the variables <ii can be 
obtained from the variables Sj by writing 



Si + 1 



(28) 



Substituting (J2HJ) into (J23]) one gets 

SgGS, V, z) = e^ v Z T {P, V, H G ), (29) 

where 
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and 
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is the Ising partition function, Eq. (|2"T|) . with 
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Now, putting into Eq. (|29[) the well-known exact for- 
mula for the partition function of the closed Ising chain 
of V spins in the external magnetic field 



\+{P,H G ) v + \-(P,H G ) v (33) 



X+W,H G ) V , 



(34) 



where 



\±(f3,H G ) = e P jG (cosh(pH G )±^smh 2 (pH G ) + e-^ J Gj, 

(35) 

the grand partition function for the one-dimensional lat- 
tice gas with nearest-neighbor interactions becomes, for 
V -> oo, 



E G (x,V,z) = {x*Jz X + (x, z)) 1 



(36) 
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where 

x = e 13 ^ 4 , z = e^, (37) 

and 

, , . X 3 \/~Z 1 X 6 Z 1 1 X 2 , 

corresponds to Eq. Q35p but written in the new variables 
x and z. Finally, inserting Eq. (|36f into Eq. (fTT|) one 
gets free energy of the one-dimensional lattice gas with 
nearest neighbor interactions in the infinite volume limit, 

<f> G (x,V,z) = -jln(xV^X+(x,z)). (39) 

Further in this paper, in Sec. IIV1 results of this section 
will be used to analyze cluster properties of the lattice 
gas model. 

IV. CLUSTERS IN LATTICE-GAS MODEL 

In this section, the combinatorial approach described 
in Sec. |TT] will be used to study properties of the one- 
dimensional lattice-gas model. In what follows, three 
cases: i. the infinite temperature, ii. the range of fi- 
nite temperatures, and iii. the zero temperature limit, 
will be discussed separately. 

A. Infinite temperature limit 

In the infinite temperature limit one has 

lim = 0. (40) 

T— too 

Therefore, since e = const one gets, see Eq. (|37)l . 

lim x = l, (41) 

T— >oo 

and the expression for the grand partition function, 
Eq. simplifies to 

lim E G (x,V,z) = {l + z) v . (42) 
Accordingly, the grand potential becomes 

lim $ G (x, V, z) = -— ln(l + z), (43) 

T->oo p 

Successive derivatives, <f> n , of the grand potential with 
respect to z and evaluated at z — 0, Eq. (J7}, give the 
following closed-form expression for the Bell polynomials' 
coefficients: 

lim w n (B) = V(n - l)!(-l) n_1 . (44) 

T— >oo 

From the last expression, it is obvious that the param- 
eters do not satisfy the conditions under which the gas 



can be considered as the ideal gas of clusters, cf. Eq. ([9])- 
By definition, lattice-gas clusters are sets of neighboring 
particles. There is no energy of interaction between such 
clusters, cf. Eq. (|2T))) . However, due to the so-called ex- 
cluded volume effect, lattice-gas clusters in some sense 
interact with each other, even if there is no direct inter- 
action between their particles. The effect has its origin 
in the cluster counting problem, and consists in the fact 
that two clusters of a given size once defined cannot ap- 
proach close enough to one another to be counted, under 
the definition, as a single cluster. Consequently, only 
at small densities, the lattice-gas can be considered as a 
perfect gas of clusters. 

Nevertheless, the case of the one-dimensional lattice 
gas in the infinite temperature limit allows to perform 
a direct validation of our combinatorial approach. By 
inserting the coefficients {u>„ (/?)}, Eq. (|44"]) . into the ex- 
pression for the canonical partition function, Eq. (|14p . 
after some algebra one gets 

N V k (-V\ N ~ k 
lim Z G (I3,V,N) = ]T — ^ B N , k ({(n-l)\}) (45) 

-t^w*»>i-G)w> 

fc=l ' v 7 

where \s(N, k) \ represents the signless (or unsigned) Stir- 
ling number of the first kind, (JQ stands for the binomial 
coefficient, and where some basic combinatorial identities 
have been used, including: i. properties of Bell polyno- 
mials ([6], pp. 133-137), i.e. 

B Ntk ({ab n x n }) = a k b N B Ntk ({x n }), (47) 

and 

B Ntk ({(n-l)\}) = \s(N,k)\, (48) 

and ii. identities involving Stirling numbers of the first 
kind Q, i.e. the relation between signed and unsigned 
Stirling numbers, 

\s(N,k)\ = (-l) N - k s(N,k), (49) 

and its generating function, 

fy\ A V k 

U=£]vtW)- ( 5 °) 

v 7 fc=l 
The obtained result, Eq. (|46|) . 

VnZoVW-Q, (51) 

is exactly the expected one. In the infinite tempera- 
ture limit, all accessible microstates are equally proba- 
ble. Therefore, the canonical partition function is just 
the number of microstates allowed, i.e. the number of 
ways to choose positions for N particles from the avail- 
able V positions. 
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B. Finite temperatures 



In the range of finite temperatures, i.e. for 



< P < oo, 



1 < X < oo, 



and for small particle densities, 

r N 1 

hm — < 1, 

y->oo V 



(52) 



(53) 



one can directly test Eq. (ITU1) , which describes the proba- 
bility that the gas of N particles, regardless of its energy, 
consists of k clusters. The expression mentioned holds 
exactly only for the ideal gas of clusters. In this case, 
however, although in general the lattice-gas does not sat- 
isfy the conditions specified by Eq. ([9]), but given small 
particle densities the coefficients w n (P) are nonnegative 
for a reasonable range of cluster sizes, n, and Eq. (fTOj) 
seems to provide good approximation for cluster statis- 
tics. 

The coefficients w n {P) of Bell polynomials in Eq. (HI 



wx{P) = V, 

w 2 {P) = V(2x 4 -3), 

w 3 (P) = V{6x s - 12x 4 
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,,(P) = V (nlx^-V - . . .) 



(54) 
(55) 
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and also the polynomials themselves can be easily cal- 
culated using Mathematica (Wolfram, Inc.). The nor- 
malized probability distribution, P(k), for the number of 
clusters, k, i.e. the right hand side of Eq. (fl"U| divided 
by Eq. (HT 



P(k) = 



B N , k ({w n (P)}) 
B N ({w n (p)}) : 



(58) 



is shown in Fig. [T] together with Monte Carlo simulations 
of the one-dimensional conserved-order-parameter Ising 
model Q , which is often used to study properties of lat- 
tice gases. 

In the figure, one can see that for N, V — const the 
average number of clusters, and consequently also the 
average cluster size strongly depend on temperature. At 
lower temperatures (i.e. for higher vales of P) particles 
try to stay in clusters. Configurations, {<7i}, with less 
number of large clusters are more likely to be observed 
with respect to configurations with a large number of 
small clusters which are typical of higher temperatures. 



C. Zero temperature limit 

1. Non-analyticity of the grand partition function 

The van Hove's theorem [H, [TcJ states that the limit- 
ing free energy, Eq. (|12p , characterizing one-dimensional 




FIG. 1: Comparison of probability distributions for the num- 
ber of clusters, k, in the one-dimensional lattice gas ob- 
tained from Monte Carlo simulations of N = 40 particles on 
V — 8000 sites (scattered points) and theoretical distributions 
P(k) given by Eq. (|58(l (solid lines). Numerical simulations 
were done using Metropolis algorithm with the coupling en- 
ergy between two neighboring particles e = 4, cf. Eq. (|20p . 
and for two different values of the inverse temperature, f3 — 1 
and 2. 



gases with short range interactions (such as the one con- 
sidered in this study), is an analytic function for all real 
positive values of T. This statement excludes phase tran- 
sitions for T > OK. The theorem, however, does not say 
anything about analyticity of the free energy at T = OK. 
In fact, the grand potential given by Eq. (|3T))) is non- 
analytic at T = OK, i.e. for 



and 



lim P = oo, 



lim x = oo. 

T-»0K 



(59) 



(60) 



Obviously, this non-analyticity is of the same mathemat- 
ical nature as the one observed in the free energy of the 
one-dimensional Ising model at zero temperature, when 
the magnetic field, H, goes to zero. 

At low temperatures, T — > OK, one finds that Eq. (|35|) 
becomes [f| 

A±09,ff G ) = e' 9jG (cosh(/3i/ G )±|sinh(/?i7 G )|(l + 0[e^ jG ])) 



_ e i3(J G ±\H G \) ^ 

Rewriting the last expression in the variables x and z, cf. 
Eq. p7|) , after some algebra one gets 



(61) 



lim X+(x, z) = xe 

T^OK 



±\ln(x 2 ^)\ 



(62) 



Correspondingly, the resulting grand partition function 
for V — > oo, Eq. (f36|). can be written as follows 



lim 

T-y0K 



Sa(x,V,z) = (>(* 2 v^)+IM*Vi)iy 



(63) 
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The expression reveals the non-analytic behavior at 

lim x 4 z = 1, (64) 

T^OK 



and 



w v {fi) = V\ x 



w 



(74) 



since 



while 



lim E G (x,V» = l, 



(65) 



lim S G (x,V,z) = x 4V z v . (66) 

C 4 Z->1 + 



From Eqs. (I?2")) - (|74")l it is evident that the considered lim- 
iting behavior meets the conditions specified in Eq. ([9]). 
It validates our cluster-based combinatorial approach to 
a general system of interacting particles and also high- 
lights its potential usefulness for the theory of phase tran- 
sitions. 



2. Cluster properties 

The non-analytic behavior of the grand partition func- 
tion just described has a very clear and convincing inter- 
pretation in terms of our combinatorial approach. Com- 
paring Eqs. (1651) and (|66)) with Eq. (JlJ one gets the fol- 
lowing expressions for the canonical partition functions: 

lim Z G (p,V,N)=5 Nf0 , (67) 

X 4 Z— > 1 ~ 

and 

lim Z G (/3,V,N) =x 4N 5 Ny , (68) 

z 4 z->l+ 

where S^j is the Kronecker delta. 

The last two expressions show that in the considered 
case of the zero temperature limit, the one-dimensional 
lattice gas with nearest neighbor interactions can be ob- 
served in only two particle configurations. The first con- 
figuration, 

Vi<j<y o-j = 0. (69) 

is consistent with Eqs. (|6"51) and (ffTTl) . and describes sys- 
tem with no particles (the so-called vacuum state). The 
second configuration, 

Vi<i<V a l = 1. (70) 

results from Eqs. (1561 and and corresponds to the 
spanning cluster of size V. 

To proceed with understanding the combinatorial ap- 
proach described in Sec. [ill it is instructive to analyze 
Eqs. (|67l) and (|6"5]l with the help of the general formula for 
the canonical partition function, Eq. (|14[) . Using proper- 
ties of Bell polynomials, i.e. 

B N<k (0,0,...,0,x j ,0,...) = S N<jk -fiftLxj, (71) 

one can show that the case when x 4 z approaches 1 from 
the left hand side, Eq. (|67)) . corresponds to: 

V„>i w n (p) = 0, (72) 

while the case when x 4 z approaches 1 from the right hand 
side, Eq. (|6"5]l . amounts to: 

V„>i ; „^y w n {0) = 0, (73) 



3. Zeros of the grand partition function 

In the zero temperature limit, i.e. for x — > oo, Eq. (|57[) 
can be approximated by 

V„>i w n ({3) = Vnlx^-V. (75) 

Inserting the last expression into Eq. (fT4")l for the canon- 
ical partition function one gets: 

1 N 

Z G (x,V,N) = — Y,B N A{Vnlx 4 ^}) (76) 

' k=l 
1 N 

= wX yk(x4){N ~ k)BNA{nl}n77) 

' fc=l 

1 N 

= jfi'E.VW-VLWk), (78) 
' fc=i 

where Eq. (|47|) has been used, and where L(N, k) rep- 
resent the so-called Lah numbers, which are defined as 
follows 

{N - 1\ N\ 
L(N,k)=( k _ 1 j—=B N , k ({ni}). (79) 

Lah numbers have an interesting meaning in combina- 
torics: they count the number of ways a set of iV elements 
can be partitioned into k nonempty linearly ordered sub- 
sets. The combinatorial meaning of L(N, k) allows a di- 
rect understanding of the canonical partition function, 
Z G (x, V, z), as given by Eq. (ff8]l. In short, Eq. (JTSJ states 
that there are L(N, k) particle configurations in which TV 
particles can be arranged into k linear clusters. All such 
configurations have the same energy, Eq. (|20"|) . and, there- 
fore, also the same Boltzmann factor, x 4 ( N ~ k \ Finally, 
the factor V k is due to cluster arrangement along the 
one-dimensional lattice, in which one assumes that ev- 
ery cluster may start at the same lattice site. Obviously, 
the arrangement factor does not take into account the 
excluded-volumc effect. For the reason it is only correct 
in the zero temperature limit. 

Now, putting the canonical partition function given by 
Eq. (|78p into the general formula for the grand partition 
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function, Eq. ((TJ), one gets 

oo TV 

E G (x,z) = l+Y / ^2- 

N=l fe=l 

Vz 



AT! 



V 



exp 



1 - a; 4 ^ 



L(N,k) (80) 



(81) 



where the infinite expansion involving Lah numbers has 
been used (see @, p. 156, or [ll[, pp. 108-113), i.e. 



exp 



tu 



1-t 



^ j>; t N u k 



iV=l k=l 



N - 1 
A: - 1 



(82) 



The grand partition function, Sq(x, z), as calculated 
above applies to the one-dimensional lattice gas in the 
thermodynamic limit. In Eq. (f5Tj) . volume plays the role 
of an extensive factor of the grand potential, Eq. (|lll) , 



/3{1 - x 4 z) 



V, 



(83) 



rather than an independent variable. It also must be 
noted that, in the considered zero temperature limit, x — > 
oo, the grand potential per unit volume, which defines 
pressure, Eq. (ITS]). 



P 



13(1- x 4 z)' 



(84) 



reveals a non-analytic dependence on z at x A z — » 1. The 
non-analyticity translates into the root of the grand par- 
tition function, Eq. (|8T]) . for z = 0. 

The significance of the zeros of the grand partition 
function was first pointed out by Yang and Lee |12j | , who 
showed that a phase transition in the sense of a non- 
analytic dependence of P on z for physical (i.e. real and 
positive) values of z can only occur when Hq(x,z) = 0. 
In the one-dimensional lattice gas with nearest neighbor 
interactions analyzed in this study, the only root of the 
grand partition function is for z — 0, i.e. it does not 
occur on the positive z axis. Therefore, one claims that 
phase transitions in the traditional sense do not exist in 
the considered gas, although its behavior in the zero tem- 
perature limit has some direct bearing to the problem of 
phase transitions. 



4- Supplementary remark 

It is interesting to see that Eq. (153")) can be derived most 
easily by inserting Eq. (|75p into Eq. (fT2j) and assuming 
that x i z < 1: 



/3x 4 



(3(1 - x 4 z) 



V. (85) 



There is a problem with proceeding in this way, however. 
At variance with our previous results, the obtained ex- 
pression is not justified for x A z > 1. There is no such 



a problem if one uses our combinatorial approach. The 
reason is that the approach described in Sec. [II] is based 
on the concept of formal power series (5, Il3j|. 

In mathematics, formal power series are a generaliza- 
tion of polynomials as formal objects. A formal power 
series is an object that just records a sequence of coeffi- 
cients. One may think of such a series as a power series 
in which one ignores questions of convergence. However, 
formal power series still allow one to employ much of the 
analytical machinery of normal power series - especially 
in settings which do not have natural notions of conver- 
gence. We believe that this perspective makes our com- 
binatorial approach to equilibrium statistical mechanics 
peculiarly well suited to handle problems with conver- 
gence, which are often encountered in the theory of phase 
transitions. 



V. SUMMARY 

The aim of this paper was to demonstrate, with a 
concrete example, how our general combinatorial ap- 
proach to interacting fluids works. To this end, the 
one-dimensional lattice gas with nearest neighbor inter- 
actions has been considered. Exploiting the mathemati- 
cal equivalence between the grand partition function for 
the gas and the canonical partition function for the one- 
dimensional Ising model in the external magnetic field, 
cluster properties of the former has been explored. Three 
cases: the infinite temperature limit, the range of finite 
temperatures, and the zero temperature limit has been 
discussed separately. In particular, in the range of finite 
temperature and for small particle densities the normal- 
ized probability distribution for the number of clusters 
has been found, and the non-analytic behavior of the 
grand potential in the zero temperature limit, having a 
direct bearing on phase transitions has been analyzed. 
The investigation of the zero temperature limit for the 
gas, has allowed us to remark on the method of formal 
power series which is behind our approach, and which, 
we believe, makes the approach peculiarly well suited to 
handle problems covered by the theory of phase transi- 
tions. 

In this study, the main purpose was to validate general 
results of our earlier paper. Apart from this purpose, 
however, we have made interesting additions to the still- 
developing theory of one-dimensional lattice gases (see 
e.g. [lj-[l7|), which has proven useful in studying many 
natural phenomena in nanophysics, surface science, and 
biophysics (see e.g. [l8l - [2l1 |). 
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